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A second order phase transition for the three dimensional 
Gross-Neveu model is established for one fermion species 
N = 1. This transition breaks a parity-like discrete symme- 
try. It constitutes its peculiar universality class with critical 
exponent v — 0.63 and scalar and fermionic anomalous dimen- 
sion rja = 0.31 and r/^, = 0.11, respectively. We also compute 
critical exponents for other N. Our results are based on exact 
renormalization group equations. 

An understanding of systems with many fermionic degrees 
of freedom is one of the big challenges in statistical physics. 
Due to the anticommuting nature of the variables numerical 
simulations are not straightforward — analytical methods are 
crucially needed. One typically has to solve a functional in- 
tegral for a d dimensional system with Grassmann variables. 
Approximate solutions for "test models" would be of great 
value. The Gross-Neveu (GN) model [1] is one of the sim- 
plest fermionic models. In three dimensions a discrete sym- 
metry forbids a mass term unless it is spontaneously broken. 
In the symmetric phase the GN model is therefore a realiza- 
tion of a statistical system of gapless fermions. For a large 
number N of fermion species it is known [2-5] that a second 
order phase transition separates the symmetric phase from an 
ordered phase where the symmetry is spontaneously broken 
and the fermions become massive. Using methods based on 
an exact renormalization group equation [6] a second order 
transition for N > 2 was confirmed. We know, however, of no 
previous work which clarifies the existence and nature of the 
phase transition in the simplest model with only one fermion 
species. The model with one fermion species is inaccessible to 
lattice simulations due to the fermion doubling problem and 
the 1/N expansion is not expected to give reasonable results 
for N = 1. The case N = 1 is also of special interest since 
an order parameter {ipj ip 3 } ^ leads to a ground state which 
does not admit any discrete symmetry involving the reflection 
of all coordinates, in contrast to the models with N > 2. 

In this letter we improve the exact renormalization group 
approach and establish a second order phase transition for 
N — 1. We also compute the critical exponents. This is im- 
portant beyond a possible relevance for real physical systems: 
the GN model constitutes a peculiar universality class due to 
the presence of massless fermions at the critical point. Just as 
the C>(iV)-Heisenberg models for bosons, the GN model could 
in the future become a benchmark for our understanding of 
critical systems in presence of fermions. 

The GN model describes N fermionic fields with local quar- 
tic interaction. Here ipj, j = 1...N, are irreducible representa- 
tions of the group O(d) including parity reflections, i. e. 2 d//2 
component Dirac spinors for d even and 2 ( - d ~ 1 ^ 2 for d odd. 
The classical Euclidean action 

S = / <l".r {«/, I r)^, • ^-(,- I r: ? ■':,)) (1) 

is symmetric under a coordinate reflection ip(x) i— > 



—ip(—x), ip(x) > ip(—x). (We note that ip and ip are in- 
dependent variables in an Eucidean formulation.) A nonva- 
nishing expectation value of ipjip 3 spontaneously breaks this 
symmetry. If the spinors ipj contain more than one irreducible 
representation of the rotation (or Lorentz) group SO(d) we 
can find alternative definitions of the coordinate reflection 
where ipjip 3 remains invariant. In particular, this is realized 
for Dirac spinors in even dimensions where a mass term cou- 
ples two irreducible representations ipL and ipR. One may 
then define a standard parity transformation under which ipip 
is invariant. In this case, however, ipRipL breaks a discrete 
chiral symmetry ipi, \— > —ipL,ipL i— ► —ipL- We will concen- 
trate here on one two-component spinor in three dimensions. 
For N — 1 the above "parity transformation" is the only 
possible choice of coordinate reflections and ipjip 3 7^ spon- 
taneously breaks this symmetry. 

Depending on the value of the coupling G the model will 
be in the symmetric phase or exhibit spontaneous symmetry 
breaking (SSB) with a nontrivial expectation value of the or- 
der parameter {Gipj ip 3 ) . We will describe a space dependent 
fermion bilinear —iGtpj{x)ip 3 {x) by a real scalar field a(x) 
such that SSB is indicated by nonzero homogeneous (a) 0. 
By a partial bosonization we express the GN model (1) as an 
equivalent Yukawa model with 

5 CT = yd^j^^^+^^+i^ 2 }. (2) 

The equivalence of the partition function can be seen 
by performing the Gaussian a-integration (rjip means 
Jd d x fjj (x) ip 3 (x) etc.) 

Z[fj, r}\ = j VaVipVip exp(— S a [o, tp, ip] + fjip + rjip) 

= JvipVipexp(^~ipi$ip + f]ip + r/ip - ^(ipip) 2 ^ (3) 
X y©aexp(-i(a + J G^) 2 ), 

where the last factor yields an irrelevant constant. 

A powerful tool for non-perturbative examinations are ex- 
act renormalization group equations for the effective action 
F, i. e. the generating functional of the 1PI Green's functions 
[7,8]. Starting with the classical action S at the UV cutoff 
A, we obtain a type of coarse-grained free energy by in- 
tegrating out the quantum fluctuations with momenta larger 
than a given scale k. Eventually, we reach the macroscopic 
thermodynamic potential [8] r at k — > 0. The IR cutoff k 
is implemented in a smooth way by introducing a mass-like 
term ASk into the classical action which gives extra masses to 
modes with momenta smaller than k. In the limit k — ► the 
IR cutoff is absent, and fluctuations at all scales have been 
taken into account. At k — > A all fluctuations are suppressed 
and Fa approaches S. 

This procedure should be explained more precisely. For no- 
tation purposes we combine the fermionic and (real) bosonic 
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fields to a column vector \ = (°'>V'>V' T )> an< i the row J con- 
tains all the external sources J — (J,fj,r) T ). We start with 
the generating functional for the connected correlation func- 
tions in the presence of the IR cutoff: 



W k [J] = In y© X cx P (S4x] - AS fc [ X ] + Jx) 



(4) 



For vanishing ASk (at k — > 0) this matches exactly with the 
free energy of the Gross-Neveu model (1). The effective action 
is then defined via a modified Legendre transformation by 



r fc [$] :=-w k [j[*]] + -AS k [$], 



(5) 



which depends on the expectation values of the fields $ = (x) • 
The infrared cutoff takes the form ASfc[$] = i$ T i? fc $, and 
Rk is a matrix: 



Rk(p,q) 




(2n) d 5 d (p- 



with R kB (q) = Z aik q 2 r k B{q) and R kF (q) = -Z^,^r kF {q). 
{Z„,k and Z^, zk are wavefunction renormalizations.) With 
these definitions, an exact renormalization group equation for 
the scale dependence of Fk can be found [7-9] : 



d t V k = isTr{(lf ' +R k r 1 <9 t R k } 



(6) 



where t = ln(fc/A). The super-trace runs over momenta and 
all internal indices and provides appropriate minus signs for 
the fermionic sector. The heart of the flow equation is the 
fluctuation matrix 



IV 



(7) 



Together with Rk it represents the exact inverse propagator 
at the scale k. 

Equation (6) is an exact but complicated functional differ- 
ential equation. There is no way around some approximation 
by truncating the most general form of T fc . We work here in 
the lowest order of a systematic derivative expansion where T fe 
contains a scalar potential, kinetic terms and a Yukawa cou- 
pling. In momentum space it is given by (fdq = f d d q/(2ir) d ): 

r fc [<7,V,$ = Jd d xU k (a) + jdq{^a(-q)q 2 a(q) 
-Z i ,^ k ipAq)i' t P J (q) + Jdpih k Tpj(p)a(p- q)^(g)j (8) 

The connection between G and h becomes clear if we rescale 
a in (2) to ha and set G = h 2 /rn 2 , m CT denoting the boson 
mass. 

Using a truncation of the effective action causes the limit 
k — > to depend on the precise form of the cutoff functions. 
In order to take control of this we have used two different 
choices, an exponential [8] and a linear [10] cutoff: 



yr% xp (y) 



y 



yr% n (y) = (l-y)Q(l-y), 



exp(y) - 1 : 

where y = q 2 /k 2 , and rp(y) is chosen in both cases such that 
2/(1 + tb) = 2/(1 + tf) 2 ■ We introduce renormalized, dimen- 
sionsless quantities h\ = Z^ 1 Z~ 2 k d ~ 4 h 2 k , p = \Z a k 2 ~ d o 2 , 

Uk = Ukk~ d , and we use u' k = ^r*- etc. 



We obtain a set of evolution equations for the effective 
parameters of the theory by inserting (8) in the exact flow 
equation (6). The evolution equation of the effective scalar 
potential u k is found by evaluating for a constant scalar 
background field: 

dtu k = ~du k + (d - 2 + r]v)pu'k 

+ 2v d {l d (u' k + 2pu'i;T hr ) -N'l[p d (2h 2 kp; V4 ,)}, (9) 

where we use N' = d-,N, with d 7 the dimension of the 7 
matrices, and v' 1 = 2 d+1 n d/2 F(d/2). The definition of the 
threshhold functions l d and l„ F ^ d as well as l n FB i2, m d ll: „ 2 , 
ra^, and m n FB J d used below (see (11)-(13)) is given in the 
appendix of [8]. These functions contain all the momentum 
integrations. For the linear cutoff the integrations can be 
performed analytically — an enormous advantage for the sub- 
sequent numerics. We obtain: 

,d, \ 2(5 n , + n)/ 7/ CT \ 1 



In WW = 



2((5„,o + n) ( ??v 



1 - 



d+1 / (I+lj)^ 1 



4^ ) 2 d (wi,W2;»7v„»?<7) = 



(1 +Wl)"l(l +UJ 2 ) ri 



1+Wi 



n»t, n2 (wi,w 2 ;77 CT ) 



+ 1 J l+cu 2 \ d+2 
1 



(1 + (x))' l l(l + 



d-2 (1 + lu) 3 \2d~4 4 J (1 + lu) 2 

(FB)d t \ I 1 V<r \ 1 

<,„> (^ 2 - M ) = " d+lj (1+ 0,0^(1 + ^ 

The anomalous dimensions jj a and 77,/, are defined as 

r? CT (fc) = -d t \nZ aik , 7?v,(fc) = -d t \nZ^ tk . (10) 

Taking second derivatives of (6) with respect to the fields 
by means of (7), we obtain evolution equations for the exact 
inverse boson and fermion propagators. Expanding the mo- 
mentum dependence at q 2 = for a — aok at the potential 
minimum yields equations for i) a and rj^ : 

r) a {k) = 8^-^K k (3u'k + 2Kku'k ) 2 mi i0 (u'k + 2n k Uk, 0; r? CT ) 

+N' h 2 k [m{ F)d (2hlK k ; 77^,) - 2h 2 k n k m 2 F) ° l (2h 2 k K k ; r/^)j j (11) 

»7V.(fc) = 8^jhlm[ F 2 B)d (2hlKk,u' k + 2nku'k;ri i ,,ri a ) (12) 

Here Kfc denotes the position of the minimum of u k (p) and 
all derivatives of the potential are evaluated at p — K k . Di- 
viding the evolution equation of the fermion propagator by a 
and evaluating it at zero momentum and a = oofc yields the 
evolution equation of h\: 

dthl = {2r]^ + r\ a + d - A)h 2 k 

+ 8v d h 4 kl[ F 1 B)d (2h 2 k Kk,u'k +2K k u'k\ri i ,,r) a ) (13) 
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We study the flow of the effective potential in two differ- 
ent ways that correspond to two different truncations of Tk- 
First, we expand the evolution equation of the potential (9) in 
a Taylor series up to third order in p around its minimum and 
evolve the resulting coupled ordinary differential equations [6] . 
In a second, more involved approach, the full equation for u' k 
is discretized on a grid, yielding coupled ordinary differential 
equations (one for each grid point) which are solved simulta- 
neously. We use the equation for u' k instead of the one for u k 
for numerical reasons. 

The only critical parameter in the theory is the four- 
fermion coupling Gact- This parameter has to be tuned in 
order to be near the phase transition. The flow equations 
(9)-(13) are evolved from an UV cutoff scale A to k — > 0. We 
concentrate on the three-dimensional case with one fermion 
species (d = 3 and N = 1). The initial values of the pa- 
rameters are chosen such that Ta = S a : Z a A = 1CF 10 ~ 0, 
Z^a = 1, h\ = cA and u' A (p) = e := (Z ct aGaA) _1 for all 
p. The GN model corresponds to c = 1 whereas for c / 1 
we investigate Yukawa type theories in the same universality 
class. 

Very close to the phase transition we find scaling solu- 
tions for all evolving parameters corresponding to vanish- 
ing beta functions. This behaviour indicates a second order 
phase transition. The symmetric regime is characterized by 
u' k (0) > as well as K k = 0, while after spontaneous symme- 
try breaking u' k (0) becomes negative and K k > 0. Starting in 
the symmetric regime the system evolves into the SSB regime 
and reaches the scaling solution. Further down the flow, it 
evolves back into the symmetric regime for Ga < Ga ci . or it 
remains in the SSB regime for Ga > Ga c ,.. This behaviour 
contrasts the one for N > 2 where the scaling solution is 
located in the symmetric regime. 




FIG. 1. Critical flow of 10n k (full line) and it' fc (0) (dashed 
line) plotted as functions of t = ln(fc/A). The thin dashed 
lines show the flow for c = 10~ 10 . We have chosen two initial 
values of c very near the phase transition. The flow separates 
only for very small k, according to the respective phase. 

Figure 1 shows the critical flow of K k and u' k (0) for two dif- 
ferent values of c, namely c = 1 (full and dashed line, resp.) 
and c = 10~ 10 (thin dashed lines). It nicely demonstrates that 
only the beginning of the flow is affected by the choice of the 



non-critical parameters. All universal quantities like critical 
exponents and mass gaps are independent of c. The poly- 
nomial expansion of the potential works well only for small 
enough c. The value of the critical coupling GAcr depends, 
however, on the choice of c. In the limit c — > 1 it converges 
towards a constant value, which is reported in table I. 

We calculate the critical exponents characterizing a second 
order phase transition. In order to test the reliability of the 
numerical algorithms, we determine the exponents in the sym- 
metric as well as in the ordered phase. The scale dependent 
renormalized boson mass rn^ R (k) is defined as the curvature 
of the potential at its minimum. In the SSB regime it is 
given by m 2 R (k) = 2k 2 K k u' k (K k ), in the symmetric regime 
by m 2 R (k) = k 2 u' k (0). In the ordered phase the running of 
m aR (k) essentially stops once k becomes much smaller than 
m,cr R (k). The situation is slightly more difficult in the sym- 
metric phase since the fermions are massless. Their fluctua- 
tions induce a scale dependence of Z a ,k even for very small 
momenta k: r) a — ► 1 for k — > in contrast to r\ a — > for 
k — > in the ordered phase. To get rid of this problem we 
define the renormalized boson mass at some fixed small ratio 
r c = k/fh aR [6] as 



- 2 

m aR 



;(n' fec (0)-u' fe c ;(0)) 



: r c m a R. 



(14) 



where u' k r (0) denotes u' k (0) on the critical trajectory. In the 
numerical calculations we have used a ratio r c — 0.01, but our 
results do not depend on r c for r c < 0.1. Thus we define the 
critical exponent v 



1 dlnrni R (Se) 
"=2^0 dlnSe ' 



(15) 



where 8e = \e — e cr \. The exponent 7 is defined as usual, since 
the unrenormalized boson mass m 2 = Z a m 2 R is not affected 
by the fluctuations of the massless fermions: 



d\nm„(8e) 
7 = hm — ^ — - 

Se^o omoe 



(16) 



In the ordered phase both exponents are defined as usual, 
using m 2 R instead of fh 2 R for the definition of v. The critical 
exponent (3 is defined as 



lim 



d In erg 



2 se^o dlnSe' 



(17) 



with (To = linifc-^o u ok . The exponents r/ a and r/^ for the criti- 
cal correlation functions are computed by taking the values of 
the scale dependent anomalous dimensions r) a (k) and J? 1 /.(fc) 
at the scaling solution (sec. 4.2 of [8]). Table I lists our re- 
sults for the critical exponents. We find a good match of the 
values in the two different phases. Besides, we have checked 
the index relations 7 = v{2 — rj a ) and (3 — \v(d — 2 + r) a ) 
which are well fulfilled. The dependence of the exponents on 
the cutoff functions rs and rF as well as on the truncation 
can also be seen from the table. The latter one seems to be 
weaker when using the linear cutoff. The error on the expo- 
nent could be larger than the difference between the results 
for different truncations of the potential and different cutoffs. 
This issue could be investigated by extending the truncation 
(8), e.g. by including a quartic fermion interaction (4>ip) 2 
along the lines discussed in [11]. 
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Since the fluctuations generate bosonic and fermionic 
masses, one might be interested in the resulting gaps. They 
are proportional to the order parameter po = \Z„oq: 

m a R = A CT p , m^ R = A^po (18) 

The gaps A CT and A^, are shown in table I. 





full eq. for u' k 


Taylor expansion 


Cutoff 


lin 


exp 


lin 


exp 


Vsymm 


0.621 


0.640 


0.623 


0.633 


^symm 


1.051 


1.077 


1.053 


1.062 


v(2-r) a ) 


1.051 


1.076 


1.054 


1.064 


Vssb 


0.620 


0.637 


0.622 


0.632 


fssb 


1.050 


1.071 


1.053 


1.062 


Pssb 


0.406 


0.420 


0.407 


0.417 




0.405 


0.420 


0.407 


0.417 




0.308 


0.319 


0.308 


0.319 




0.112 


0.114 


0.112 


0.113 


A(j 


16.0 


17.6 


16.8 


18.1 


Ay, 


14.2 


14.9 


14.4 


15.2 


GrAcr A 


43.13 


26.68 







TABLE I. Critical exponents and mass gaps for N = 1 in 
three dimensions 



We have also calculated the critical exponents in the three- 
dimensional GN model for N > 1. They are compared with 
results from other methods in table II. For N = 2, 4 and 12 
critical exponents have been calculated in the 1/A-expansion 
to 0(1/N 2 ) [12,13] and anomalous dimensions to 0(1/N 3 ) 
using conformal techniques [14]. Monte-Carlo methods have 
been used to calculate critical exponents for = 4 which are 
compared with results from the e = 4 — d expansion to 0(e 2 ) 
[13]. Our values for the Taylor expansion with exponential 
cutoff (b) agree well with [6] where the precise numerical im- 
plementation was different. 

The largest discrepancy between different approaches con- 
cerns the values of the anomalous dimensions. This is similar 
as for the bosonic O(N) models and is expected to improve 
substantially in the next order in the derivative expansion 
[15]. The overall picture is, however, quite satisfactory and 
lends support to the validity of our method based on the exact 
renormalization group equation. We believe that our finding 
of a second order phase transition for N = 1 is quite robust. 



IN 


Z 


A 

4 


1 




0.927™ 


1.018" 


1.018™ 


V 


0.962 6 


1.016 6 


1.011 6 




0.738 c 


0.903 c 


1.007 c 






0.948™ 








1.00(4) e 






0.525™ 


0.756 


0.927 


n<? 


0.554 6 


0.786 6 


0.935 6 




0.635 c 


0.776 c 


0.914 c 






0.763 


0.913-f 






0.754(8) 






0.071 


U.032 


0.0087 




0.067 b 


0.028 i> 


0.0057 6 


0.105-f 


0.044^ 


0.0124 J 


A CT 


15.1° 


12.0" 


5.5 




17.3 6 


14.0 6 


6.4" 


Ay, 


10.8° 


7.5" 


3.3" 




11.4 6 


7.9 6 


3.4 6 



TABLE II. Critical exponents for various N in three di- 
mensions 

(a) linear cutoff, full equation for u' k 

(b) exponential cutoff, Taylor expansion of u' k 

(c) 1/A-cxpansion (Pade-Borel resummed) 

(d) e-expansion (Pade-Borel resummed) 

(e) Monte-Carlo simulations 

(f) conformal techniques 
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